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ABSTRACT 

A new wavefront sensing approach, derived from the successful curvature wavefront sensing 
concept but using a non-linear phase retrieval wavefront reconstruction scheme, is described. 
The non-linear curvature wavefront sensor (nlCWFS) approaches the theoretical sensitivity limit 
imposed by fundamental physics by taking full advantage of wavefront spatial coherence in the 
pupil plane. Interference speckles formed by natural starlight encode wavefront aberrations with 
the sensitivity set by the telescope's diffraction limit A/D rather than the seeing limit of more 
conventional linear WFSs. Closed- loop adaptive optics simulations show that with a nlCWFS, 
a 100 nm RMS wavefront error can be reached on a 8-m telescope on a my = 13 natural 
guide star. The nlCWFS technique is best suited for high precision adaptive optics on bright 
natural guide stars. It is therefore an attractive technique to consider for direct imaging of 
exoplanets and disks around nearby stars, where achieved performance is set by wavefront control 
accuracy, and exquisite control of low order aberrations is essential for high contrast coronagraphic 
imaging. Performance gains derived from simulations are shown, and approaches for high speed 
reconstruction algorithms are briefly discussed. 

Subject headings: instrumentation: adaptive optics — techniques: high angular resolution 



1. Introduction 

Adaptive Optics (AO) systems are now an es- 
sential part of ground-based telescopes, and rou- 
tinely deliver diffraction limited images at near- 
IR and mid-IR wavelengths on 8 to 10 m tele- 
scopes. For several key science applications, the 
required wavefront quality is however still well be- 
yond what current systems deliver. For example, 
obtaining high quality diffraction-limited images 
at visible wavelengths requires residual wavefront 
errors to be below 100 nm RMS. Direct imaging 
of exoplanets or circumstellar debris disks with a 
high contrast coronagraphic camera is even more 
demanding, and requires exquisite control of low 
order aberrations and mid-spatial frequencies. 

In the last decade, several high performance 
coronagraphs concepts have been developed in 
the laboratory, and full AO+coronagraph lab- 



oratory systems have demonstrated that both 
deformable mirror and coronagraph technolo- 
gies are now sufficiently mature to allow direct 
imaging of planets down to Earth mass provided 
that wavefront aberrations are small: for exam- 
ple, the vacuum High Contrast Imaging Imaging 
Testbed (HCIT) a NASA JPL has achieved better 
than 10~ 9 raw PSF contra st at 4 A/D separa- 
tion (iTrauger fc Traubl 120071 ) . This high contrast 
value is orders of magnitude better than, at a few 
X/D, the w 10 -3 best raw PSF contrast currently 
achieved in the near-IR on 8 m to 10 m tele- 
scopes and the expected « 10 -4 raw PSF contrast 
for Extreme- AO syste ms currently under assem- 
bly such as SPHERE (jBeuzit et al.ll2008l) on the 
Very La rge Telescope (VLT) a nd Gemini Planet 
Imager ((Macintosh et alj [2008) . The large gap 
between laboratory demonstrations and on-sky 
performance is entirely due to wavefront control 
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accuracy. A key difference between laboratory 
systems and on-sky conditions is the timescale 
of wavefront variations, which requires fast wave- 
front sampling and therefore fundamentally limits 
the number of photons available for sensing on 
ground-based systems. 

The wavefront information which can be ex- 
tracted from a given number of photons is known 
from first principles, and accurately defines the 
"ideal" performance which can be reac hed by 
AO s ystems on ground-based telescopes (jGuvon 
2005) . Direct comparison between this ideal limit 
and the sensitivity offered by current wavefront 
sensing schemes reveals a huge performance gap. 
Commonly used wavefront sensors such as Shack- 
Hartmann and Curvature are very robust and flex- 
ible, but are poorly suited to high quality wave- 
front measurement: for low-order aberrations on 
8 to 10 m telescopes, they require ~ 100 to 1000 
times more photons than more optimized wave- 
front sensing techniques. While they offer nearly 
ideal sensitivity at the spatial frequency defined 
by the subaperture spacing, they suffer from poor 
sensitivity at low spatial frequencies, which are 
most critical for direct imaging of exoplanets and 
disks. 

This paper describes a high sensitivity wave- 
front sensing approach which is optically derived 
from the curvature wavefront sensing technique, 
and which shares the same non-linear reconstruc- 
tion algorithm as used for phase diversity. In 
phase diversity, images near the focal plane are 
used for wavefront reconstruction with a non- 
linear algorithm, while curvature wavefront sens- 
ing uses images near the pupil plane and a linear 
wavefront reconstruction algorithm. Phase diver- 
sity is theoretically more sensitive than linea r cur- 
vature wavefront sensing ( Fienup et al.lll998 ). but, 
for ground-based adaptive optics applications, suf- 
fers in practice from smearing of focal plane speck- 
les due to chromatic effects and time averaging of 
the fast-changing PSF. I propose in this paper a 
solution to improve curvature WFS sensitivity by 
using non-linear signals in the near-pupil images, 
and offer a solution to optically mitigate the chro- 
matic smearing of the speckles which would other- 
wise reduce sensitivity. Unlike phase diversity, the 
proposed scheme is robust against time averaging 
of wavefront aberrations, which gradually brings 
the proposed technique closer to a less sensitive 



but stable linear curvature wavefront sensor. 

The non-linear curvature wavefront sensor (nl- 
CWFS) principle is described in EJSJ The wavefront 
reconstruction algorithm for the nlCWFS is de- 
scribed in fj3l The WFS performance is evaluated 
through numerical simulations in 2] and discussed 



2. nlCWFS Principle 

Figure [1] shows the proposed optical imple- 
mentation for the non-linear Curvature WFS (nl- 
CWFS). In this section, the main differences from 
a conventional curvature WFS are explained: 

• The nlCWFS achieves its sensitivity by us- 
ing non-linear signal in images that are ac- 
quired in planes which are considerably fur- 
ther away from the pupil, where non-linear 
effects are strong ( §2.ip 

• While curvature systems use a variable fo- 
cal length vibrating membrane, here, beam 
splitters and lenses are used to produce the 
required conjugation planes on the detec- 
tor (s) array f ^2~2]) 

• Chromatic lenses are used in the sensor's fo- 
cal planes for compensation of chromatic ef- 
fects ( g23 ) 

• Four conjugation planes are created instead 
of two in conventional curvature ( fl2.4p 

2.1. Non-linear signal and large defocus 
distances 

The performance of conventional Curvature 
WFS (cCWFS) is strongly constrained by the re- 
quirement that wavefront sensor signals must be a 
linear function of th e input wavefront aberrations 
( Guyon et al. 2008f ). This linear constraint keeps 
the defocus distance much smaller than it should 
be for optimal conversion of phase aberrations 
into intensity fluctuati ons. The dua l strok e lin- 
ear scheme proposed in lGuvon et al. ] (l2008h only 
slightly mitigates the negative effects of these con- 
straints, but linearity is still enforced. 

A far more serious limitation is that linear re- 
construction fully ignores the non-linear coupling 
between high and low order aberration in the de- 
focused pupil images. This is best visualized by 
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Fig. 1. — Conceptual layout of a non- linear Curvature wavefront sensor. See text for details. 



direct comparison of linear vs. non-linear mea- 
surement of tip-tilt from defocused pupil images. 

In cCWFS, tip-tilt is measured as an overall 
displacement of the pupil in the defocused im- 
ages: only the edges of the defocused pupil im- 
ages are used to extract signal. The displace- 
ment is proportional to the defocus distance dz. 
If the wavefront is otherwise perfectly flat, tip-tilt 
is theoretically best sensed by an extremely large 
(infini te) defocusing distance ( van Dam fc Lang 
l2002bh . where the defocused pupil image is in 
fact close to a focal plane image, and the tip-tilt 
measurement sensitivity is equivalent to what is 
achieved when measuring tip-tilt from a properly 
sampled X/D- wide focal plane diffraction- limited 
image. This is hardly realistic, as higher order 
wavefront aberrations also need to be sensed, and 
the defocus distance must therefore be kept small, 
and not all photons should be allocated to tip- 
tilt sensing. Even if the defocus distance were 
very large, wavefront aberrations would preclude 
diffraction-limited imaging at the wavefront sens- 
ing wavelength. The defocused pupil plane images 
are therefore convolved by the A/ro-wide atmo- 
spheric seeing kernel (were tq is the Fried param- 
eter quantifying the strength of the atmospheric 
turbulence): unless the AO correction is quite 
good, the tip-tilt sensitivity is at best equivalent to 
measuring the centroid of a A/ro-wide spot. It is in 
fact not as good as that, as the defocus distance is 
finite, and tip-tilt is therefore only measured from 
detectors at the edge of the pupil. 

If non-linear effects are to be taken into account 
for tip-tilt measurement, the situation is radically 
different. First, larger defocus distances are al- 



lowed, which would improve the tip-tilt measure- 
ment sensitivity even if it were derived with the 
simple "pupil photocenter" diagnostic (which is 
linear). The real benefit however comes from the 
fact that the defocused pupil images are "speck- 
led" due to edge (pupil edge, central obstruction, 
spiders) diffraction effects and high order aberra- 
tions. These highly contrasted speckles, visible 
in figure [2l are X/D wide (which means they are 
physically smaller closer in to the pupil plane, and 
infinitely small - invisible - at the pupil plane). 
Provided that the defocused pupil image is ac- 
quired with proper spatial sampling, a tip-tilt can 
therefore be measured as the displacement of a 
"cloud" of A/D-wide speckles. High order aberra- 
tions must be known to some degree to be able to 
recover this highly sensitive tip-tilt measurement. 
Similarly, low order aberrations need to be known 
to measure high order aberrations from the defo- 
cused pupil images. 

2.2. Beam splitters 

In conventional curvature systems, the required 
defocus is achieved sequentially with a resonant 
vibrating membrane (VM) in the system's focal 
plane. Because the high performance of the pro- 
posed wavefront sensing technique relies on re- 
solving speckles in the defocused pupil images, 
care must be given to avoid smearing them, which 
would happen with a sine-wave motion of the 
membrane. If a VM were used, its motion should 
therefore be "squared", and the short transition 
period during which the defocus altitude is rapidly 
changing should not be used toward wavefront 
sensing (this time could conveniently be available 
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Fig. 2. — Simulated frames obtained by the nlCWFS in monochromatic light (top), in polychromatic light 
with no chromatic compensation (center) and in polychromatic light with chromatic compensation (bottom) . 
The wavefront error across the 8m diameter pupil used for this simulation is shown on the right. 



for detector readout). 

A better solution is to simultaneous image the 
four defocused pupil images on a single detector 
array (or on four smaller detector arrays) with a 
small number of beam splitters and mirrors. The 
desired beam size on the WFS detector array is 
on the order of 1mm (100 pixels accross the pupil 
for 10 /xm pixels). In this scaled down copy of 
the telescope beam, a 1000 km vertical propaga- 
tion distance is scaled down to approximately 1 
cm. In the non-linear curvature wavefront sens- 
ing scheme described in this paper, there is also 
no need to adapt the defocus altitudes as a func- 
tion of atmoshperic conditions, so such a "static" 
implementation appears suitable. 

2.3. Chromatic compensation with the fo- 
cal plane lenses 

The primary purpose of lenses L1+, L1-, L2+ 
and L2— in the WFS's focal planes is to achieve 
the desired altitude conjugation on the detector 
planes: if the lenses are removed, four identical 
pupil plane images would be acquired. In the nl- 



CWF, these lenses arc designed to be chromatic 
so that, in each of the 4 channels, the conjugation 
altitude is inversely proportional to wavelength. 
Figure [2] illustrates this chromatic compensation. 
Without compensation, speckles in the polychro- 
matic images are slightly blurred. This is more 
apparent in with larger propagation distances, for 
which the finest-sized speckles disappear. With 
the chromatic compensation, polychromatic defo- 
cus pupil images are very close to monochromatic 
images, and the blurring of the smallest speckles 
is very minimal. For very large propagation dis- 
tances, defocused pupil images are close to a focal 
plane image, and the proposed chromatic compen- 
sation is equivalent compensating for the A scaling 
of the position of PSF diffraction features in the 
focal plane. 

2.4. Number of defocus planes 

In lGuvon et al. I ( 2008h . the dual stroke sensing 



scheme was justified by the need to mitigate some 
of the constraints brought by the requirement of 
preserving a linear relationship between WFS sig- 
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nals (curvature) and input wavefront phase. With 
a non-linear wavefront reconstruction scheme, it 
would therefore seem that the need for dual stroke 
is removed. 

With the non-linear reconstruction algorithm 
and simulation parameters presented in this pa- 
per, the AO loop could not be closed using two 
symmetric detector planes. No attempt has been 
made to explore other configurations (for exam- 
ple, 2 non symmetric planes, or 3 planes) or sta- 
bilize the loop by modifying or constraining the 
reconstruction algorithm, so this does not exclude 
the possibility that nlCWFS may work with fewer 
than 4 planes. 

While the reasons for the loop instability in a 
symmetric 2 planes nlCWFS have not been stud- 
ied, it may be due to a fundamental property of 
diffractive propagation. The nlCWFS achieves its 
sensitivity by acquiring images at large propaga- 
tion distances z from the pupil plane. In such an 
image, there is a spatial frequency / = y/2/{zX) 
for which z is the Talbot propagation distance at 
which the intensity modulation vanishes. It may 
therefore be important to measure the pupil in- 
tensity at two different values of z to ensure that 
every spatial frequency within the control range of 
the AO system produces an intensity signal in the 
WFS. 

2.5. Detector sampling 

The effect of detector sampling for nlCWFS has 
not been numerically quantified in this work, but 
physical considerations suggest it is driven by two 
requirements: 

1. The non- linear dual stroke curvature achieves 
its high sensitivity for low order aberrations 
by "tracking" X/d wide speckles in the defo- 
cused pupil images. The sampling therefore 
needs to be sufficiently fine to image these 
speckles with 2 pixels per X/D for maxi- 
mum sensitivity. The linear number of pix- 
els across the pupil should therefore be at 
least N = 2D 2 /(Xz) where z is the defocus 
distance and D is the telescope diameter. 
With D = 8m, X = 0.8^m and z = 2000km, 
we find that the pupil should be 80 pixels 
across. Coarser sampling can be achieved by 
pushing the defocus altitudes higher, with 
little or no loss in performance since the re- 



construction algorithm can accurately take 
into account non-linear effects. 

2. As the defocus distance is increased, Fresnel 
diffraction "blurrs" the pupil image, which 
becomes larger than the geometric pupil. 
This effect is visible in fig. [2] This is es- 
pecially important for the highest spatial 
frequencies which the WFS must measure, 
since high spatial frequencies will diffract 
at larger distances away from the geomet- 
ric pupil. The "blurred" pupil size becomes 
Dbiurred = D + 2za where a is the greater 
of two values: a\, the angular "control ra- 
dius" of the AO system at the sensing wave- 
length, or a.2, half the angular size of the 
PSF at the wavefront sensing wavelength. In 
closed loop, (Xi is usually much smaller than 
the natural seeing, and ot2 is usually smaller 
than at. With D = 8m, z = 2000km, 
A = 0.8/im, and a = 0.2"= 10 X/d (ap- 
proximately 300 modes sensed), the detector 
needs to extend almost 2 m in every direc- 
tion around the edge of the 8 m geometric 
pupil size. 

When accounting for both effects, the linear num- 
ber of pixels across the detector should be: 

N= 2D{D + 2za) _ 
Xz 

This equation clearly shows that smaller values of 
z require more pixels. To minimize the number 
of pixels read, it is best to chose z equal to or 
greater than zq = D/(2a), which is the propaga- 
tion distance for which the "blurred" pupil size is 
twice the geometric pupil size. For D = 8m and 
a = 0.2", z = 4000km. For z = z , N = N = 
4Da/X, and for z » z , N = = N /2. 

With z >> zq, the total number of pix- 
els required per conjugation plane is N^ir/4 = 
irD 2 a 2 / X 2 , which is equal to the number of X/D- 
wide speckles within the control radius a of the 
AO system. Since each of these focal plane speckes 
corresponds to two degrees of freedom (amplitude 
and phase), a pair of defocused pupil images con- 
tains as many pixels as there are modes sensed by 
the nlCWFS. With two such pairs of defocused 
images read, and if z > Zq, the nlCWFS therefore 
requires at least 2 pixels to be read per WFS mode. 
Compared to a Shack-Hartmann WFS with 2x2 
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pixels cells or a pyramid WFS designed to sense 
the same total number of modes, the non-linear 
dual stroke curvature WFS is therefore expected 
to require as many or slightly more pixels to be 
read. 

Equation [T] has not been verified by numerical 
simulations, and the non- linear iterative nature of 
the reconstruction may produce a more compli- 
cated picture than the one presented in this sec- 
tion. Both the required detector sampling and the 
required number of conjugation planes are impor- 
tant drivers for system cost and complexity, and 
will need to be explored in the future by compar- 
ing closed loop numerical simulations for several 
combinations of defocus distances, pixel scale for 
each plane, and detector size for each plane. Re- 
construction algorithms should also be chosen to 
mitigate coarse detector sampling or limited num- 
ber of conjugation planes if possible. 

3. Wavefront reconstruction algorithm 

In this section, an iterative algorithm which 
successfully reconstructs the wavefront from the 
four defocused pupil images is identified. While 
more efficient algorithms may exist for this prob- 
lem, this algorithm is then used in the rest of this 
paper to estimate the nlCWFS performance. 

3.1. Phase retrieval 

The wavefront may be reconstructed from sev- 
eral defocused pu pil images through a phase di- 



versity algorithm (Gonsalvcs & Chidlaw 



1979h . A 



Gershberg-Saxon ( Gerchberg fc Saxtonl 197 j ) al- 



gorithm is chosen in this work for simplicity and 
flexibility, but many other opti ons and variat ions 
on this algorithm are available (jFienup 1119821 ), as 
much work has been devoted to phase retrieval 
algorithms in the last few decades. Similar phase 
retrieval algorithms have been successfully used to 
measure o ptical aberrations in the Hubble Space 
Telescope ilFienup et alll993tlRoddier fc Roddieil 
1993t iKrist fc Burrowslll995h . The algorithm used 
for this work is illustrated in figure [H 

Using this algorithm, nearly optimal closed 
loop performance has been obtained in simulations 
with as few as 10 GS loops by (1) fast forwarding 
the convergence thanks to the fact that in con- 
secutive GS iterations, the solution moves in the 
same direction, and (2) running a few lower res- 



olution GS iterations on binned data to rapidly 
approach the solution prior the full size slower 
GS iterations. Essentially all of the computation 
time is taken by the full-size FFTs in the propa- 
gations, while the phase unwrapping described in 
the next section is very quick, as it is both lin- 
ear and non-iterative (see §3.2|l . Execution speed 
is a steep function of sampling of the defocused 
pupil images. Eight complex 64x64 pixels 2D 
FFTs are required per GS iteration (assuming a 
45 pixel diameter pupil). On a single modern 
CPU, each FFT takes 0.088 ms (number adopted 



from http://www.fftw.org/speed/ for a 3 GHz In 



tel Core Duo CPU), corresponding to 7ms total 
FFT time. On a single CPU, the iterative nl- 
CWFS discussed in this section can therefore run 
at «100 Hz with current hardware. 

Other approaches to this non-linear prob- 
lem e xist, but have not been tested in this 



work. Ivan Dam fc Land (|2002q ) have suggested 



a geometry-based algorithm to derive wavefront 
aberrations from defocused image. The recon- 
struction can also be linearized in a w lrad do- 
main around any chosen wavefront, and it should 
therefore be possible to effectively combine a 
coarse non-linear algorithm with a fine resolution 
linear algorithm. 

3.2. Phase unwrapping and filtering 

The output of the GS iterative algorithm is 
a phase map of the pupil, which then needs to 
be unwrapped and filtered before commands are 
sent to the DM. Unwrapping may not be neces- 
sary if the closed loop performance is such that 
the peak-to-valley phase is less than 2ir. Figure 
[3] shows how unwrapping and filtering (in this ex- 
ample, projection on a set of Zernike polynomial) 
can be done together in one step. The phase slope 
is insensitive to phase wrapping effects, as it is 
here computed as difference between phase of ad- 
jacent pixels, and any value above 7r or below — 7r 
is brought back in the — 7r to n interval by adding 
or subtraction 2tt. For each of the two axis (x 
and y) the resulting phase slope is decomposed 
as a linear sum of precomputed phase slopes (one 
per Zernike polynomials). The coefficients of this 
decomposition are the Zernike coefficients of the 
wavefront, which can therefore simultaneously be 
unwrapped and filtered by being reconstructed as 
a sum of Zernike polynomials. If the output of 
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Fig. 3. — Proposed wavefront reconstruction algorithm for non-linear curvature wavefront sensing. 



the GS algorithm is sufficiently good, the differ- 
ence between the wrapped pupil phase and the 
filtered/unwrapped pupil phase (this difference is 
small enough to be free of phase wrapping effects) 
may be added to the final unwrapped phase map 
to get back high spatial frequencies removed by 
the filtering step. If however the WFS frames are 
very noisy, some filtering may also be included 
within the GS loops to force the GS algorithm to 
produce a smooth phase map, which will mitigate 
problems associated with phase unwrapping errors 
on noisy wavefront data. The final step of the 
wavefront estimation is to appropriately weigh the 
output of our wavefront reconstruction algorithm 
with prior knowledge of the wavefront and knowl- 
edge of the statistical properties of the wave front 
(|Gendron fc Lenalll994t llrwan fc Lanelll998f ) and 
its temporal evolution. This optimization prob- 
lem has not been explored in this paper, and is 
not specific to the nlCWFS wavefront sensor (al- 
though the exact "tuning" of this optimization is 
wavefront sensor-dependent). 



4. Performance: Numerical Simulations 

In this section, numerical simulations are per- 
formed to verify and evaluate the nlCWFS 's per- 
formance in a closed loop adaptive optics system. 
In fc|4.H reconstruction of a single wavefront is 
tested using the algorithm shown in Fig. [3l Closed 
loop simulations of nlCWFS and SH based AO sys- 
tems are presented in §4.2l to illustrate the sensitiv- 
ity gain offered by nlCWFS. Chromaticity effects 
are investigated in ^4.31 with and without the op- 
tical chromaticity compensation proposed in this 
work. In g3 the flux limit at which nlCWFS 
starts to loose the ability to see individual speck- 
les is probed through numerical simulations. Fi- 
nally, in ^4.51 closed loop simulations are shown 
for sparse apertures. 

4.1. Reconstruction of a single aberrated 
wavefront 

Figure [4] illustrates the performance of non- 
linear CWFS. For this simple simulation, the ini- 
tial wavefront aberrations were 609nm RMS across 
the pupil. This level of aberration is, for a 8m 
telescope, intermediate between what is expected 
to be encountered in open loop (typically 2000nm 
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Fig. 4. — Wavefront reconstruction using the algorithm shown in fig. [3) Four noisy defocused pupil images 
(images (a), (b), (c) and (d)) are acquired to transform the pupil phase aberrations (e) into intensity signals. 
The input pupil phase is 609 nm RMS, yielding the PSF (f) before correction. After correction, the residual 
pupil phase aberration (g) is 34.4 nm RMS, allowing high Strehl ratio imaging (h). All images in this figure 
were obtained at 0.65 /im. The total number of photons available for wavefront sensing in 2e4. 



RMS) and in closed-loop (around lOOnm RMS for 
a bright guide star) . The number of photons avail- 
able for wavefront sensing was chosen to be 2e4, 
which corresponds to approximately 0.5 ms inte- 
gration time on a my = 10 source in a 0.5/im wide 
bandwidth and a 20% system efficiency (prod- 
uct of optics throughput and detector quantum 
efficiency). This example is therefore very rep- 
resentative of situations commonly encountered 
in ground-based Adaptive Optics when observing 
moderately bright targets. 

Results shown in figure [4] demonstrate the nl- 
CWFS's good sensitivity: with only 2e4 photons 
(photon noise is very pronounced in the WFS 
frames) , the wavefront is reconstructed to a 34nm 
RMS accuracy. The PSF obtained is therefore 
diffraction limited in the visible. In the next part 
of this paper, simulations are performed in condi- 
tions closer to a real AO system: in §4.31 the effect 
of chromaticity is evaluated, and in §4.21 and 34.41 
detailed closed-loop simulations are used to verify 
performance. 



4.2. Closed-loop simulation: performance 
and direct comparison with SHWFS 

Figure [5] compares closed loop AO performance 
obtained with SHWFS and nlCWFS for a total 
flux of 3xl0 6 ph/s on a 8-m telescope (magnitude 
13 with 20% efficiency and a 0.5 ^m wide band- 
pass). Parameters adopted for this simulation are 
given in Table Q] Both the WFS and imaging 
wavelength are 0.85 /im. The WFS detector array 
is considered perfect (no readout noise, no dark 
current). For the SH systems, each subaperture 
PSF is imaged on 64x64 pixels, and the DM is 
considered ideal, correcting all spatial frequencies 
measured. Four SHWFS configurations were sim- 
ulated with D/d varying from 60 to 9 (D = tele- 
scope diameter, d = subaperture size). For the 
nlCWFS simulation, the detector acquires simul- 
taneously the four defocused pupil images with 90 
pixels across the pupil in each image. The DM is 
assumed to perfectly correct all modes up to 16 
cycles per aperture (CPA), and provides no cor- 
rection for higher spatial frequencies. 
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Table 1 

Default parameters used for simulations 



nlCWFS SHWFS 



Telescope diameter 8m 

Seeing 0.6"FWHM at 0.5 /mi 

Wind speed 10 m.s~ 

WFS wavelength 0.85 /im 

WFS number of subapertures not applicable from 9 to 60 subapertures across pupil 

WFS defocus distances -3500 km, -2500 km, 2500 km, 3500 km not applicable 

WFS detector sampling 90 pixels across pupil 64x64 pix. per subaper. (pix. size <C A/d) 

WFS detector readout noise 

AO loop gain 1.0 

AO loop modal filtering none (gain = 1 for all modes within control range) 

Spatial frequency control range up to 16 cycles per aperture up to CPA = D/3d 

Edge subapertures not applicable active if > 70 % illuminated 



The simulations code was written in C lan- 



guage, and used the fftw library (jFrigo k, Johnson 



20051 ) for fast Fourier transforms. In each case, AO 
loop frequency was optimized to minimize resid- 
ual wavefront aberration (RMS phase error over 
the pupil): this optimization was done by run- 
ning the loop at different frequencies in 20Hz in- 
crements and selecting the optimal solution. Fig- 
ure [5] illustrates the effect of noise propagation in 
a SHWFS: with a large number of subapertures 
(D/d = 60), sensing accuracy for low order aber- 
rations is poor, and the overall wavefront RMS 
error is large even though a high number of modes 
is corrected. With poor wavefront sensing sensi- 
tivity, the optimal loop frequency is low (140 Hz 
for D/d=60) to balance time lag and photon noise. 
With fewer larger subapertures (D/d = 36), cor- 
rection is better for low spatial frequencies, and 
residual wavefront error drops even though high 
spatial frequencies are not corrected, as can be 
seen from the light just outside the darker square 
indicating the control region of the AO system. 
Further reducing the number of subaperture how- 
ever fails to significantly improve correction of low 
spatial frequencies as the SHWFS sensitivity is 
now dominated by seeing size (the subapertures 
are larger than tq). For a given guide star bright- 
ness, there is an optimal number (higher number 
for brighter stars) of subapertures in a SHWFS for 
which sensitivity loss on low order aberrations is 
balanced against number of modes measured. 

Thanks the nlCWFS 's high sensitivity across 
the full range of spatial frequencies measured, the 



PSF halo is much fainter than could be achieved 
with any of the SHWFS configurations tested. 
This improved sensitivity also allows the nlCWFS- 
based AO system to run faster (260 Hz) . Increas- 
ing the number of wavefront modes measured does 
not lower performance for low order aberrations in 
a nlCWFS. 

4.3. Chromaticity 



As described in §2.3[ the proposed design nl- 
CWFS uses chromatic lenses to improve the wave- 
front sensor's ability to work in broadband light. 
This scheme does compensate for chromaticity in- 
duced by diffraction propagation, but will blur the 
defocused pupil images in the geometrical optics 
regime where aberrations are large. For example, 
a large amount of tip/tilt a will shift defocused 
pupil images by adz independently of wavelength. 
The chromatic lenses proposed in §2.31 will blur 
the images since dz is now a function of A. It is 
therefore important to verify that the use of chro- 
matic lenses does not prevent the loop from ini- 
tially closing when aberrations are large, and that 
the correction remains stable in time. 

Simulation results shown in Figure [H] confirm 
that the chromatic lenses scheme allows the nl- 
CWFS to both close the loop and maintain a 
good correction in broadband light with little loss 
of performance compared to the monochromatic 
case. Simulations shown in Figure [5] use the same 
parameters as used in the nlCWFS result shown 
in Figure [5j except for dzl, which was changed 
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Fig. 5.— 
nlCWFS. 



Simulated 10 second exposures closed loop PSFs delivered by AO systems with SHWFS and 



from 2500 km to 300 km. Reducing dzl does help 
with loop stability in polychromatic light, at very 
little cost in performance. Without the chromatic 
lenses, the loop would not have been stable in a 
20% wide band (Figure bottom). 

These results are very encouraging, as they were 
obtained by using a wavefront reconstruction al- 
gorithm which assumes the nlCWFS is working 
in monochromatic light. It therefore appears that 
polychromatic operation of a nlCWFS requires no 
additional computing power. 

Broadband use of the nlCWFS also requires an 
atmospheric dispersion compensator (ADC) ahead 
of the nlCWFS. For nlCWFS to operate at full ef- 
ficiency, chromaticity introduced by atmospheric 
dispersion should be accurately compensated to 
keep the lateral PSF shift below the diffraction 
limit of the telescope across the wavefront sens- 
ing spectral coverage: even a small angular extent 
of the source used for wavefront sensing is highly 
detrimental. 



4.4. Coherence loss and optimal Wave- 
front Sensing wavelength 

Closed loop simulations using the nlCWFS 
show a rapid drop in performance at about m «15 
(Figure[8|). This limit is due to a loss of coher- 
ence in the pupil plane: the number of photons 
becomes too low to track the diffraction-limited 
speckles which are essential to the nlCWFSs sen- 
sitivity. Beyond this limit, the nlCWFS operates 
in a similar regime as conventional linear WFSs, 
and delivers partially corrected PSFs (Figure [7|). 

The results obtained in figure O also show that 
for faint guide stars (approximately 14 < my < 17 
in the figure), it is preferable to choose a longer 
wavefront sensing wavelength, thanks to the re- 
covery of coherence in the pupil plane. For bright 
stars, a shorter sensing wavelength is better, since 
the same OPD will correspond to a larger phase 
shift, which can be detected more easily. Figure 
[8] does not show this last point because the RMS 
residual wavefront error for bright stars is domi- 
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Fig. 8. — Simulated performance of a low order nlCWFS-based system (first 200 Zernike modes corrected) 
as a function of sensing wavelength (0.7, 0.85 and 1.0 /im) and guide star brightness. The stellar magnitudes 
given in this figure assume a 20% efficiency in a 0.5/xm wide band. Each point in this figure is a 1-second 
average of the RMS residual wavefront error. See text for details. 



nated by the number of modes corrected in this 
"low-order" AO system. 

4.5. Sparse apertures, spiders and low or- 
der aberrations 

For both existing and future telescopes, the 
pupil is usually not a single continuous disk as 
assumed so far in this work. Narrow spider vanes 
have a very small effect on most wavefront sensing 
schemes, but thicker spider vanes can be a more 
serious problem, especially if a Shack-Hartmann 
sensor is used with subaperture size comparable 
to or smaller than the spider thickness. This sit- 
uation will likely occur on the next generation of 
large telescopes. This problem is even more se- 
rious for a large telescope composed a few large 
segments, with large gaps between the segments, 



of which the Giant Magellan Telescope is the per- 
fect example. Incorporating statistical knowledge 
of the phase error can help, but will not work 
for large gaps or vibration-induced piston between 
segments. A separate scheme to measure differen- 
tial piston between segments is often required to 
solve this problem. 

To illustrate this problem, I now consider a tele- 
scope with 20cm thick spider vanes, and a Shack- 
Hartmann wavefront sensor. To improve the AO 
system performance for bright stars, a large num- 
ber of small subapertures is preferred. With lm 
subapertures, the spiders would not be a serious 
problem, but with 20cm subapertures, the sub- 
apertures along the spiders are either fully dark 
or partially illuminated. No subaperture sees light 
simultaneously from both sides of the spider and 
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Fig. 9. — Wavefront reconstruction with dual stroke non-linear CWFS on a sparse pupil. The pupil amplitude 
(a) and phase (b) yield the four defocused pupil images shown on the right. The recovered pupil phase (c) 
and the residual phase error (d) demonstrate that dual stroke non-linear CWFS can simultaneously measure 
OPD within and across segments. This polychromatic simulation was performed with 2e8 photons in a 
d\/\ = 0.5 wide band centered at 0.65 fim. 



the WFS is therefore blind to an OPD between the 
two sides of the spider. Even if the subapertures 
are twice the spider thickness, the amount of light 
which is mixed between the two sides of the gap is 
relatively small, and OPD measurement between 
the two sides will be noisy. The ability to mea- 
sure wavefront across a gap is a function of how 
much light is mixed within the WFS between the 
two sides of the gap. Shack-Hartmann sensors, 
especially with small subapertures, are therefore 
poorly suited for this task. 

As illustrated in figure M with a nlCWFS the 
images acquired are sufficiently far from the pupil 
plane to introduce significant mixing between 
pupil "segments", therefore enabling high sensi- 
tivity measurement of OPD between segments. A 
fundamental problem however persists if the gaps 
are wide: no monochromatic wavefront sensing 
scheme can detect OPDs which are multiples of 
A. The solution to this problem (which has not 
yet been explored) is to intelligently use chro- 
matic effects: a trial and error iterative algorithm 
using the focal plane science image (which, pre- 
sumably is at a different wavelength) may be used 
to "lock" the inter-gap OPD within A - once this 
is achieved, the closed loop nlCWFS-based AO 



system will accurately maintain the correct inter- 
gap OPD. Alternatively, the broadband defocused 
pupil images acquired within the WFS should be 
sufficient to solve for this problem: the "inter- 
ference" fringes between the segments are highly 
contrasted if the inter-gap OPD is close to zero, 
but much less contrasted if it is close to A. Us- 
ing this as a diagnostic could reveal the absolute 
value of the intergap OPD, not its sign: some it- 
erative trial-and-error algorithm seems necessary 
to initially lock the AO control loop in the correct 
intergap OPD. 

The fundamental gain (over a Shack-Hartmann 
sensor) offered by non-linear CWFS for sensing 
inter-gap OPD measurement is merely a demon- 
stration of its higher sensitivity for low-order aber- 
rations sensing. An inter-gap OPD is mostly a 
low-order aberration. In both cases (inter-gap 
OPD and low order aberration), the key to high 
sensitivity lies in being able to mix light from rel- 
atively distant parts of the pupil. 
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Fig. 10. — Optical sensitivity of several wavefront sensor options. The number of photon (y axis, relative) 
required to achieve a fixed wavefront sensing SNR is shown as a function of spatial frequency (x-axis). 
Sensitivities are shown for wavefront quality ranging from perfectly flat (top left) to 20 rad RMS (bottom 
right) . 
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Fig. 6. — Top: Closed loop 1 second expo- 
sure PSFs at 0.85 fim with a nlCWFS work- 
ing in monochromatic and broadband (dA/A = 
40%) light. In both simulations, the total num- 
ber of photons in the nlCWFS is 3e6 ph/s. Bot- 
tom: Without chromatic correction within the nl- 
CWFS (replacing the chromatic lenses by achro- 
matic lenses in Figure [1} , the loop would be un- 
stable with dA/A=0.2. 

5. Discussion 

5.1. Sensitivity comparison with Other 
Adaptive Optics Wavefront Sensing 
Techniques 

Figure [10] shows for different wavefront sens- 
ing schemes how strong an optical signal is pro- 
duced when small changes are introduced in the 
wavefront of a 8 m diameter telescope. The WFS 
is operating at 0.85 /im, and all phases are mea- 
sured at this wavelength. The "strength" of the 
signal is measured as signal to noise (SNR) for 
a fixed number of available photons, or, equiva- 
lently, how many photons are required to reach 
a fixed SNR. Although only photon noise is con- 
sidered, the highest sensitivity WFSs will also be 
more robust against detector readout noise. The 
WFS schemes compared in Figure fTOl are: 



Fig. 7. — Simulated long exposure 1.6 ^m PSFs 
obtained with a non-linear Dual stroke Curvature- 
based AO system. The sensing wavelength is 0.85 
/xm for this simulation. 

• Two SHWFSs with different number of sub- 
apertures (D/d — 16 and D/d = 32), both 
using linear wavefront reconstruction. 

• A curvature WFS with adz = lOOOfcm prop- 
agation distance on either side of the pupil 
plane and a linear wavefront reconstruction 



A lin ear dual-stroke curvature WFS (jGuvon et al. 
2008) with dzl = 800fcm and dz2 = 2000fcm 



• A linear pyramid wavefront sensor where the 
pyramid position is fixed (Fixed Pyramid 
WFS, or FPYRWFS). 

• A FPYRWFS where non-linear signal is 
taken into account for wavefront reconstruc- 
tion sensitivity. 

• A curvature WFS where non-linear signal is 
taken into account, with two defocus dis- 
tances dzl — 2500/cm and dz2 — 3500fcm 
(4 pupil images total). 

The y-axis shows how many photons are needed to 
achieve a fixed wavefront sensing signal-to-noisc 
ratio (SNR) compared to the o ptimal sensitiy - 
ity wavefront sensor as defined in iGuvonl (|2005h . 
For example, the top left panel (0 rad RMS case) 
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shows that at 1 cycle per aperture (CPA) , both the 
SHWFS with D/d = 32 and the curvature WFS 
require 200 times more photons than an ideal 
WFS to produce an intensity signal with the same 
SNR. This flux factor is given as a function of spa- 
tial frequency (x axis), measured here in CPA. Fig- 
ure [10] was generated by (1) Computing the inten- 
sity response in the WFS, (2) Computing how this 
intensity response changes when a small aberra- 
tion is added to the wavefront, and (3) Comparing 
this change, pixel by pixel, to the photon noise and 
combining the resulting SNRs into a single over- 
all SNR. This computation was done not only for 
a flat incoming wavefront, but also for aberrated 
wavefronts, and therefore quantifies optical sensi- 
tivity beyond the small pert urbation linear regime 
considered in lGuvon (2005). A distinction is also 
made between linear schemes (SH, conventional 
curvature, conventional pyramid), where the in- 
tensity signals from a series of uncorrelated wave- 
fronts were added, and non-linear schemes, where 
a single aberrated wavefront is used. This compu- 
tation step is necessary to average out non-linear 
signals in linear WFSs: for example, a SHWFS 
cannot use the shape and speckles of individual 
spots to unambiguously measure high order aber- 
rations within each subaperture; and the centroids 
measurement error is driven by the overall size 
of the spots instead of the size of speckles within 
these spots 

Figure[Tn]shows why a Curvature WFS in which 
non-linear signal is used appears to be the most 
promising approach. Other existing WFS designs 
cannot as easily be modified to offer similar bene- 
fits: 

• The SHWFS suffers from poor sensitivity at 
low spatial frequencies. This poor sensitiv- 
ity is fundamentally due to the splitting of 
the pupil into subapertures which prevents 
interferences between distant points in the 
pupil. 

• The linear fixed pyramid WFS is very sen- 
sitive at all spatial frequencies in the small 
aberration regime, but shows poor sensitiv- 
ity at low spatial frequencies when the wave- 
front quality is less than ideal. This orig- 
inates from both the linear constraint and 
an optical range limitation due to the fact 
that the inner parts of the PSF can spend 



a significant fraction of the time completely 
in a single quadrant of the pyramid - in this 
configuration, the WFS does not produce a 
strong optical signal for low order aberra- 
tions. If the non-linear part of the signal 
is used, sensitivity is improved, but is still, 
at low spatial frequencies, about an order of 
magnitude poorer than it would be for a flat 
wavefront. 

• The linear curvature WFS lacks sensitivity 
at low spatial frequencies. The dual-stroke 
scheme offers a modest (factor 2) sensitivity 
improvement at low spatial frequencies when 
wavefront erro rs are moderate. As shown in 

pupil 



Guvon et al. 



images, mixing of light between distant parts 
of the pupil occurs, but it requires large de- 
focus distances to be easily seen and it is 
highly non-linear. If this non-linear signal is 
taken into account to derive the WFS sensi- 
tivity, the dual-stroke CWFS is highly sensi- 
tive at all spatial frequencies when the wave- 
front is highly distorted. Curiously, simula- 
tions show that sensitivity is poorer when 
the wavefront quality is good: in this case, 
the speckles which carry the non-linear in- 
formation in the defocused pupil images are 
weak or absent. In this regime, a static 
known aberration should be added in the 
WFS to maintain full sensitivity. A 2 A RMS 
static aberration with a Kolmogorov power 
spectral density has been included in Figure 
[TUlin the rad and 4 rad cases. 

This analysis therefore shows that the non- 
linear curvature wavefront sensor (nlCWFS) is, 
among all the approaches listed, the only one to 
offers near-ideal optical sensitivity at all spatial 
frequencies and all levels of wavefront aberrations 
within the to 20 rad RMS range considered. The 
simulations shown in this section only considered 
the effect of a single wavefront mode at a time, and 
the results obtained do not automatically mean 
that the nlCWFS (or non-linear pyramid) signals 
can unambiguously be used to accurately measure 
the wavefront. For example, different wavefront 
aberrations could produce identical nlCWFS sig- 
nals. The existence of a wavefront reconstruction 
algorithm which can recover a complex wavefront 
map was shown in iJJl 
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5.2. Comparison with phase diversity 

In Phase Diversity (PD), images near the fo- 
cal plane are acquired to measure wavefront aber- 
rations. Similarly to nlCWFS, a non-linear re- 
construction algorithm is used to estimate pupil 
phase from several images acquired in different 
planes. The position of these planes is usually 
given in pupil defocus phase for PD while it is 
given in propagation distance for nlCWFS. Physi- 
cally an image of the pupil plane acquired after a 
free space propagation is equivalent to a defocused 
focal plane image. The two quantities are related 
by the following equation: 



Tld 2 

AXl 



(2) 



where dcj> is the pupil defocus term in radian peak- 
to-valley d is the telescope diameter and / is the 
propagation distance. 

In conventional linear curvature wavefront sens- 
ing, the propagation distance I is about 1000 km 
or less, equivalent to at least 10 waves (peak to 
valley) of defocus for PD. PD operates at smaller 
defocus values, typically 1 to 5 waves, and it has 
been shown that operating PD with larger defocu s 
values reduces its sensitivity ([Fienup et al.lll9981 ). 

In the nlCWFS example adopted for Figure 
[TOl I is 3500 km, yielding d<fr = 2.3 waves (peak 
to valley). nlCWFS therefore acquires images 
which are physically very similar to PD images, 
and erases much of the fundamental difference 
previously identifie d between PD and curvature 
wavefront sensing ( Fienup et al. 19981 ). The re- 



sults presented in this paper, and especially Fig- 
ure [T01 are consistent with the sensitivity compar- 
ison bet ween PD and c urvat ure previoulsy pub- 
lished in iFienup et al. ( 19981). One key finding 
of their study is that curvature WFS's sensitiv- 
ity is fundamentally limited by the large defocus 
value (= small propagation distance in the pupil 
pl ane) at which it ope rates (as shown in Figure 5 



Fienup et alj (j.1998)). This limitation is in fact 



due to the linearity constraint in conventional cur- 
vature, but is lifted in the nlCWFS scheme which 
allows non-linear signals to b e used at large prop - 
agation distances. Figure 5 of lFienup et al. (19981) 
shows that optimal sensitivity is reached when the 
defocus is below w 5 waves for a wavefront aber- 
ration including the 52 lowest order Zernike poly- 
nomials. With higher order aberrations included 



and a well sampled detector, the same figure would 
show that the defocus can be increased further 
(the propagation distance can be reduced) because 
the acquired images would still contain the X/D- 
wide speckles necessary for high sensitivity wave- 
front sensing. 

5.3. nlCWFS performance for Extreme- 
AO on bright targets 

5.3.1. nlCWFS linearity in small aberration 
regime 

When wavefront variations are small (less than 
1 rad), the defocused pupil images acquired by 
the nlCWFS become less contrasted and they are 
dominated by diffraction ringing from the pupil's 
edges: a known static aberration should then be 
added to the beam in the WFS path to maintain 
full sensitivity (see §5.1|) . If the wavefront aberra- 
tions around this static wavefront offset are small 
(<1 rad), the nlCWFS is also linear. 

In this low aberration regime, the nlCWFS is 
then ideally suited for Extreme- AO: 

• A very fast linear wavefront reconstruction 
can be employed 

• The nlCWFS is achromatic and can work in 
broadband 

• The nlCWFS can measure with full sensitiv- 
ity wavefronts which are far from flat. The 
nlCWFS can be operated in a linear control 
loop around a "reference" wavefront which 
is not flat. This feature is very useful when 
non-common path errors need to be removed 
for example. 

Thanks to the nlCWFS linearity computing re- 
quirements can be greatly relaxed for ExAO: a 
single matrix multiplication would be performed 
to reconstruct the wavefront. The linearity rela- 
tionship also facilitates spreading the computing 
load between CPUs, which can be done more eas- 
ily than for the iterative algorithm proposed in 
Figure [3] 

5.3.2. Expected performance in a coronagraphic 
Extreme- AO system 

The performance of a nlCWFS-based extreme- 
AO system c an be quantifi ed using the same for- 
malism as in iGuvonl (|2005h . In this section, only 
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the PSF halo contributions due to time lag in the 
AO loop and photo n noise i n the WFS are con- 
sidered (term C 2 in iGuvonl (|2005h ). Scintillation 
and optical pathlength chromaticity between sens- 
ing and imaging wavelength are not taken into ac- 
count. 

Figure QTJ shows, for a photon rate of 7.8 10 7 
ph.s -1 in the WFS (my w 7 star observed in 
a 0.1 urn. wide band with no loss in the optics 
and an ideal detector), how the PSF halo surface 
brightness C*2 varies with both spatial frequency 
and WFS sampling time. With a SHWFS with 
D/d = 50, PSF halo surface brightness is domi- 
nated by photon noise in the WFS for sampling 
rates beyond 1 kHz (Figure ITT)) : increasing the 
WFS sampling rate beyond 1 kHz offers no bene- 
fit. The nlCWFS's high sensitivity allows perfor- 
mance (lower PSF halo surface brightness) with 
higher WFS sampling rates. The performance dif- 
ference is especially large at small angular separa- 
tion, where a full order of magnitude can be gained 
if the WFS can be sufficiently fast. 

6. Conclusions 

Simulations and analysis presented in this pa- 
per show that the non linear Curvature WFS tech- 
nique is potentially able to work close to the fun- 
damental theoretical sensitivity limit imposed by 
photon noise even if the PSF is not diffraction lim- 
ited at the WFS wavelength. It is also a highly 
flexible WFS choice, as it can efficiently sense 
wavefront on sparse pupils, and performs an ex- 
plicit wavefront reconstruction, therefore allow- 
ing open-loop operation and compensation of non- 
common path errors. 

In the preliminary study presented in this pa- 
per, encouraging results were obtained on interme- 
diate brightness guide stars with closed loop simu- 
lations using a relatively simple control algorithm. 
Further performance gains are to be expected if 
both modal control (optimal weighting of wave- 
front modes according to their measurement SNR 
and variance in the atmospheric turbulence) and 
predictive control are also included. A key chal- 
lenge to fully take advantage of this technique is 
to achieve sufficiently high AO loop frequency, and 
algorithms with faster convergence than the itera- 
tive loop used in this paper should be developed. 

The nlCWFS appears especially well suited for 



Extreme-AO systems aimed at delivering nearly 
flat wavefronts to a high contrast imaging cam- 
era (coronagraph) for direct imaging of exoplanets 
and disks. In these systems a nlCWFS can de- 
liver for low order aberrations on a 8m telescope 
the same wavefront measurement accuracy as a 
SHWFS with w 1000 times fewer photons. The 
nlCWFS can take full advantage of the telescope's 
diffraction limit for wavefront sensing, with a sen- 
sitivity scaling as D 4 (combination of a D 2 gain 
due smaller diffraction limit and a D 2 gain due to 
the telescope's collecting power) while more con- 
ventional WFSs scale only as D 2 (telescope's col- 
lecting power). This difference is especially large 
for the next generation of Extremely Large Tele- 
scopes (ELTs), where the use of high sensitivity 
WFS techniques such as the nlCWFS may allow 
direct imaging of exoplanets in reflected light. 

The author is thankful to Roger Angel, Phil 
Hinz, Michael Hart, Markus Kasper, Marcos van 
Dam and Christophe Verinaud for valuable com- 
ments on this manuscript. 
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